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ABSTRACT 



Aims. The goal of this work is to investigate the power spectrum of the magnetic field associated with the giant radio halo in the 
galaxy cluster A665. 

Methods. For this, we present new deep Very Large Array total intensity and polarization observations at 1.4 GHz. We simulated 
Gaussian random three-dimensional turbulent magnetic field models to reproduce the observed radio halo emission. By comparing 
observed and synthetic radio halo images we constrained the strength and structure of the intracluster magnetic field. We assumed that 
the magnetic field power spectrum is a power law with a Kolmogorov index and we imposed a local equipartition of energy density 
between relativistic particles and field. 

Results. Under these assumptions, we find that the radio halo emission in A665 is consistent with a central magnetic field strength 
of about 1.3 yuG. To explain the azimuthally averaged radio brightness profile, the magnetic field energy density should decrease 
following the thermal gas density, leading to an averaged magnetic field strength over the central 1 Mpc 3 of about 0.75 /jG. From the 
observed brightness fluctuations of the radio halo, we infer that the outer scale of the magnetic field power spectrum is ~ 450 kpc, and 
the corresponding magnetic field auto-correlation length is ~100 kpc. 
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In the hierarchical model of structure formation, clusters of 
galaxies undergo several merger events during their lifetimes. 
Shocks and the turbulence associated with a major cluster 
merger event are thought to accelerate particles and compress 
magnetic field in the intracluster medium (e.g. Roettiger et al. 
1999). 

The existence of relativistic particles and a magnetic field 
is expected to lead to large-scale diffuse emission, which is not 
associated with individual host galaxies but is a general prop- 
erty of the intracluster medium. Radio observations reveal the 
presence of sources known as halos. Radio halos are faint (~ 
lyi/Jy arcsec~ 2 at 1.4 GHz), steep-spectrurrQ (a >1) synchrotron 
sources located at the center of galaxy clusters (e.g. Feretti & 
Giovannini 2008, Ferrari et al. 2008). As recently shown by 
Cassano et al. (2007) and Murgia et al. (2009), these diffuse ra- 
dio sources can have quite different length scales, but the largest 
halos are the most luminous, in such a way that all these sources 
may have similar synchrotron emissivity. 

About 30 radio halos nowadays are known, and all of them 
are found to be strictly related to intense merger activity (e.g. 
Schuecker et al. 2001; Buote 2001; Govoni et al. 2001). Indeed, 
the most luminous radio halos are mainly associated to hot and 
massive clusters. Cluster merger events are expected to release 
a significant amount of energy in the intracluster medium. This 
energy is injected on large spatial scales, and then turbulent cas- 
cades may be generated (e.g. Brunetti & Lazarian 2007). 
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In the past years, progress has been made by analyzing the 
state of turbulence of the intracluster medium. Signature of tur- 
bulence have been found by analyzing the gas pressure of the 
Coma cluster, where the pressure fluctuations are consistent 
with a Kolmogorov power spectrum (Schuecker et al. 2004). 
Turbulence in the intracluster medium has also been studied in 
the radio band through Faraday rotation measures. Recent re- 
sults show that the magnetic field power spectrum can be esti- 
mated if very detailed rotation-measure images of cluster radio 
galaxies are available (EnBlin & Vogt 2003; Murgia et al. 2004; 
Vogt & EnBlin 2005; Govoni et al. 2006; Guidetti et al. 2008; 
Laing et al. 2008). Typically, rotation measure images of cluster 
radio galaxies permit investigating the fluctuations of the intra- 
cluster magnetic field below a spatial scale of about 50 4- 100 
kpc. On the other hand, radio halo images reveal that the intra- 
cluster magnetic field is spread over Mpc scales. Thus, it would 
be important to investigate the magnetic field turbulence over 
such a large volume of space. In fact, although radio halos are 
typically unpolarized, in A2255 (Govoni et al. 2005) and MACS 
J0717. 5+3745 (Bonafede et al. 2009), a polarized signal associ- 
ated to the radio halo has also been detected. These observations 
indicate that, at least in these clusters, the intracluster magnetic 
fields can be ordered on scales of a hundred kpc. 

In this paper we present a study of the magnetic field power 
spectrum in the galaxy cluster A665, which contains a known 
radio halo. Following Tribble (1991) and Murgia et al. (2004), 
we propose to study the cluster magnetic field strength and struc- 
ture on the basis of the radio halo properties. Indeed, informa- 
tion on the cluster magnetic field can be derived from detailed 
images of the radio halo, since the halo brightness fluctuations 



2 V. Vacca et al.: The intracluster magnetic field power spectrum in Abell 665 

Table 1. Details of the VLA observations of Abell 665. 



Obs. pointing, RA 


Obs. pointing, DEC 


Obs. frequency, v 


Bandwidth 


Config. 


Date 


Time 


Project 


(J2000) 


(J2000) 


(MHz) 


(MHz) 






(h) 




08:30:53.000 


+65:51:30.000 


1465, 1415 


25 


C 


24-M-2005 


13 


AG690 


08:30:54.830 


+65:49:54.800 


1365, 1435 


50 


D 


14-M-1996 


5 


AF304 



and the polarization level are strictly related to the intracluster 
magnetic-field power spectrum. For example, lack of polariza- 
tion and a smooth and regular surface brightness may indicate 
that the cluster magnetic field is tangled on smaller scales than 
the resolution of the radio images, while a disturbed radio mor- 
phology and the presence of polarization could be related to a 
magnetic field ordered on larger scales than the observing beam 
(Govoni et al. 2005). 

We investigate the total intensity and the polarization proper- 
ties of the A665 radio halo by means of a new deep VLA obser- 
vation at 1 .4 GHz in C configuration and previously published 
data at 1.4 GHz in D array (Giovannini & Feretti 2000). The 
power spectrum of the intracluster magnetic field fluctuations is 
constrained by comparing the data with the expectations of sim- 
ulated random 3-dimensional magnetic field models character- 
ized by different power spectra. 

The structure of this paper is as follows. In Sect. 2 we give 
a brief description of the cluster of galaxies A665. In Sect. 3 we 
present the radio observations and describe the data reduction. 
In Sect. 4 we develop the magnetic field modeling. In Sect. 5 we 
report the simulation results and the radio halo analysis. In Sect. 
6 we analyze the depolarization of the discrete radio sources. 
Finally, in Sect. 7 we summarize our conclusions. 

Throughout this paper we adopted a ACDM cosmology with 
H = 71 km &- 1 Mpc-\ Q. m = 0.27, and Q A = 0.73. At the dis- 
tance of A665 ( z=0.1829, Gomez et al. 2000), 1" corresponds 
to3.04kpc. 



2. The cluster of galaxies Abell 665 

The cluster of galaxies A665 is well known to contain an ex- 
tended diffuse radio halo first observed by Moffet & Birkinshaw 
(1989) and successively confirmed by Jones & Saunders (1996). 
At an angular resolution of about 1', the radio halo is asymmet- 
ric because it is elongated in the southeast-northwest direction 
(Giovannini & Feretti 2000). 

The total halo flux density at 1 .4 GHz is 43 . 1 + 0.8 mJy, while 
at 0.327 GHz the flux density is 197 + 6mJy. This results in an 
averaged spectral index of ahi = 1.04+0.02 (Feretti et al. 2004). 
Like all clusters containing a radio halo, A665 shows signatures 
of a major merger. Optical (e.g. Gomez et al. 2000) and X-ray 
(e.g. Hughes & Birkinshaw 1994) studies have revealed a dis- 
turbed dynamical state. 

Chandra satellite data have provided a detailed gas tempera- 
ture map revealing a shock region with a temperature jump from 
8 keV to 15 keV, located near the southern boundary of the radio 
halo (Markevitch & Vikhlinin 2001; Govoni et al. 2004). The 
radio morphology seems to follow the X-ray elongation, trail- 
ing the shock wave direction. A spectral analysis of Chandra 
data by Million & Allen (2009) indicates the possibility of non- 
thermal (or quasi-thermal) X-ray emission, likely due to inverse 
Compton scatter between relativistic electrons responsible for 
radio halo emission and cosmic microwave background photons. 



3. Radio observations and data reduction 

We present a new deep, polarization sensitive, full-synthesis ob- 
servation of A665 at 1.4 GHz obtained using the Very Large 
Array (VLA) in C configuration. The details of the observation 
are summarized in Tablefl] The data reduced following standard 
procedures using the NRAO's Astronomical Image Processing 
System (AIPS) package. 

The nearby phase calibrator 0834+555 was observed over 
a wide range in parallactic angle to separate the source polar- 
ization properties from the feed parameters. The radio source 
3C286 was used both as primary flux density calibrator and 
as reference for the absolute polarization angles. Radio inter- 
ferences were carefully excised, and several cycles of self- 
calibration and CLEAN were applied to remove residual phase 
variations. A circular beam of 15" was used to restore the fi- 
nal images of total intensity / and Stokes parameters Q and U. 
Images of polarized intensity P - ^Q 2 + U 2 (corrected for the 
positive bias) and fractional polarization FPOL — P/I were de- 
rived from the /, Q, and U images. The / image has an rms noise 
level of 07 25 /Jy/beam, while the P image has an rms noise 
level of crp 15 /iJy/beam. 

In order to improve Mv-coverage and sensitivity to the dif- 
fuse emission but still keeping a good angular resolution, we 
combined the new C configuration data with the D configuration 
data by Giovannini & Feretti (2000), see TableQ] The D con- 
figuration data were been calibrated in phase and amplitude by 
using the sources 0841+708 and 3C48, respectively. A circular 
beam of 55" was applied to restore the final total intensity im- 
age /, which has a noise level of 60 yuJy/beam. The combined 
C+D data were then self-calibrated to produce a final total in- 
tensity image. This image was restored with a circular beam of 
25" and has a noise level of 45 yt/Jy/beam. Since the D configu- 
ration data were not polarization sensitive, no Q and U images 
were produced with the combined C+D data set. The total inten- 
sity images at different resolution of the radio emission in A665 
are shown in Fig.Q] 

The central regions of the galaxy cluster are permeated by a 
low-surface brightness radio halo which is surrounded by many 
discrete radio sources. The full extent of the radio halo is readily 
visible in the D configuration image shown in Fig.QJ. However, 
the comparatively low resolution of this image does not per- 
mit easily distinguishing the contribution of the discrete radio 
sources from the diffuse radio halo. 

The new C configuration image at a resolution of 15" is 
shown in Fig.QJ). Here we can clearly separate the halo emis- 
sion from the discrete radio sources. In particular, there are a 
few embedded sources in the central regions, and at least part of 
the northwestern elongation of the halo seen at lower resolution 
comes from a blending of several discrete sources distributed, in 
projection, along an arc-like structure. At the sensitivity level 
of our C configuration observation, we did not detect polar- 
ized emission from the radio halo. Since the innermost regions 
of the halo have an average brightness of / 100 //Jy/beam, 
we can derive a lcr upper limit to the fractional polarization of 
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Fig. 1. (a): Total intensity radio contours of A665 at 1.4 GHz (VLA in D configuration) with an FWHM of 55" x 55". The first 
contour level is drawn at 180 yuJy/beam, and the rest are spaced by a factor V2. The sensitivity (lcr) is 60 /iJy/beam. (b): Total 
intensity radio contours of A665 at 1.4 GHz (VLA in C configuration) with a FWHM of 15" x 15". The first contour level is drawn 
at 75 /iJy/beam, and the rest are spaced by a factor V2. The sensitivity (lcr) is 25 yuJy/beam. (c): Total intensity radio contours of 
A665 at 1.4 GHz (combining VLA data in C and D configuration) with an FWHM of 25" X 25". The first contour level is drawn at 
135 /dy/beam, and the rest are spaced by a factor V2. The sensitivity (lcr) is 45 /Jy/beam. The contours of the radio intensity are 
overlaid on the Chandra X-ray image. The adaptively smoothed X-ray image is in the 0.8-4 keV band (Govoni et al. 2004). 



FPOL < Jcr 2 p /(I 2 - cr]) ^ 15%. The polarization of the dis- 
crete sources in the field of view is discussed in Sect|6] 

Finally, In FigOJ we present the C+D configuration image 
at a resolution of 25"(corresponding to 75 kpc at the cluster dis- 
tance). The radio iso-contours^ are overlaid on the adaptively 



smoothed Chandra X-ray image in the 0.8-4 keV band (Govoni 
et al. 2004). The X-ray and radio peaks do not overlap precisely. 
The distance between them is about 45". The radio halo emis- 
sion peaks at RA^OS'^l''^ 1 and DEC=65 fl '50 m 2F, while the X- 
ray peak coordinates are RA=08 /l 30 m 59 5 and DEC=65 fl '50 m 21 s . 



2 The reader can refer to Fig.f4^ for an alternative visualization with 
contours and color scale of the radio emission at 1.4 GHz. 



4 



V. Vacca et al.: The intracluster magnetic field power spectrum in Abell 665 



4. Magnetic field modeling 

The halo brightness fluctuations and the polarization level can 
be related to the intracluster magnetic-field power spectrum. 
Following the approaches proposed by Tribble (1991) and 
Murgia et al. (2004), we analyzed the radio images presented in 
the previous section to constrain the intracluster magnetic field 
in A665. We proceeded in two steps. First, we performed a 2- 
dimensional analysis to determine the radial scaling of the mag- 
netic field, i.e. the average magnetic field strength at the cluster 
center and its radial decline. Second, we focused on the radio- 
halo total intensity and polarization fluctuations with respect to 
the average radial profiles to determine their power spectrum 
with the aid of 3-dimensional numerical simulations. 



4. 1 . The radial scaling of the equipartition magnetic field 
strength 

In the following 2-dimensional analysis, we want to infer the av- 
erage magnetic field strength at the cluster center and possibly 
its scaling with the thermal gas density (see also Murgia et al. 
2009). The observed radio halo intensity is directly related to 
the integral of the magnetic field and relativistic particle distri- 
bution along the line-of-sight. However, since the synchrotron 
emissivity traces the product of electron and magnetic field en- 
ergy densities, disentangling the two contributions is not possi- 
ble from radio observations alone. In this work we assume that 
the electron energy spectrum is a power law of the form 



N(e, 6) = N e~ 



, sin 6 



(1) 



with the energy e ranging from y,„i„m e c 2 and 7 max m e c 2 , with 
J max » Jmin- Such an energy spectrum is expected in the steady 
state of a continuous injection of power-law distributed cosmic- 
ray electrons subject to radiative losses. The index b is related to 
the observed spectral index throughout the relation 6 — 2a + 1 . 
The relativistic electrons are supposed to have an isotropic distri- 
bution of the pitch angle, 8, between the direction of their veloc- 
ity and the magnetic field. We set 6 = 3, in agreement with the 
observed spectral index ahi - 1 (see Sect.|2]i. We set the normal- 
ization No supposing a perfect equipartition condition between 
magnetic field and relativistic particle energy density at each 
point in the intracluster medium. In this case the radio source 
is on a minimum energy conditional and the radio emissivity at 
1 .4 GHz is given by 

J lACHz - 8.12 x 10- 43 (y m , ;l /100)^ G (ergs-'cm^Hz- 1 ). (2) 

In Eq.|2]we took the isotropy of the electron population into 
account and we averaged over all the possible directions between 
the magnetic field and the line-of-sight, i.e. the field is supposed 
to be completely tangled on an infinitesimally small scale (for a 
more realistic model see Sect. 14. 2\ . We furthermore assumed that 
the magnetic field strength scales as a function of the thermal gas 
density according to 



B(r) = Bo 



'n e (r) 
v « e (0) 



(3) 



3 The energy densities of the magnetic field, u B , and cosmic ray par- 
ticles, u cr , are exactly equal if the minimum energy condition is calcu- 
lated by assuming a fixed energy range for the cosmic ray and (5 = 3, 
as early noted by Brunetti et al. (1997) and as we did in this work and 
in Murgia et al. (2009). If, instead, fixed frequency range is used one 
recovers the classical results that u B = (3/4)u cr (Pacholczyk 1970). 



where Bq is the average magnetic field strength at the cluster cen- 
ter, while n e (r) is the thermal electron gas density profile taken 
to follow the /3-model (Cavaliere & Fusco-Femiano 1976): 



n e (r) = n e (0) \ 1 + ^ 



2\-T 



(4) 



where n e (0) is the thermal gas density at the cluster center and 
r c is the cluster core radius. In the case of A665, the gas density 
parameters have been taken as derived by Roussel et al. (2000) 
from ROSAT X-ray data and rescaled to our chosen cosmology 
( B = 0.763, r c = 112" and n c (0)=3.25- 10~ 3 cirT 3 ). 

Under these hypotheses, the synchrotron emissivity scales as 



J\ACHz(r) = Jo 1 + 



-m 



(5) 



Since radio halos are optically thin sources, the radio brightness 
profile results from the integral along the line-of-sight of the syn- 
chrotron emissivity in Eq.|5] which gives 



/(r x ) = / (l + -^ 



-3/3)7+0.5 



(6) 



where I(/ ± ) is the brightness value at the r ± projected distance 
from the cluster center, and Io is the central radio halo brightness. 
By fitting the model in Eq|6]to the data, we obtain the value of 
Io and 77, from which Bo can be estimated. We calculated the 
equipartition magnetic field by neglecting the energy contribu- 
tion from relativistic protons and other, non-emitting, heavy par- 
ticles, i.e. we used k — (where k is conventionally the ratio of 
the heavy particle to electron energy densities). Adopting k > 
leads to a higher equipartition magnetic field strength: for 6 = 3, 
B oc (1 + a:) - 25 , see also Beck & Krause (2005). In this respect, 
the magnetic field values computed in this way should be con- 
sidered as lower limits. 

In Figf2]we show the azimuthally averaged radio-halo bright- 
ness profile obtained from the C+D configuration image at 
25"resolution. Each data point represents the average brightness 
in concentric annuli of half beam width centered on the X-ray 
peak, as shown in the inset. Discrete sources have been masked 
out and excluded from the statistics. The observed brightness 
profile is traced down to a level of 3cr 7 and the best fit of the an- 
alytical model in Eq.|6]is showed. The fit was performed in the 
image plane as described in Murgia et al. (2009). To properly 
take the resolution into account, the model in Eq.|6]was first cal- 
culated in a 2-dimensional image, with the same pixel size and 
field of view as that observed, and then convolved with the same 
beam by means of a fast Fourier transform. The resulting image 
is masked exactly in the same regions as for the observations. 
Finally, the model is azimuthally averaged with the same set of 
annuli as were used to obtain the observed radial profile. All 
these functions were performed at each step during the fit proce- 
dure. As a result, the values of the central brightness, Io, and the 
index 77 provided by the fit are deconvolved quantities, and their 
estimate includes all the uncertainties related to the masked re- 
gions and to the sampling of the radial profile in annuli of finite 
width. 

The best fit of the analytical model yields a central bright- 
ness of /o = 0. 76+0 ^ yuJy/arcsec 2 and 77 = 0.47+ ' )3 , the corre- 
sponding equipartition magnetic strength at the cluster center is 
B = 1.58+°-[£ yuG. Indeed, the value of 77 we found is very close 



to the physical situation in which the magnetic field energy den- 
sity scales as the thermal gas density: B 2 oc n e . 
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Fig. 2. Analytical fit to the observed radio-halo brightness ra- 
dial profile. Each point represents the azimuthally averaged ra- 
dio brightness at 1.4 GHz obtained in annuli, as shown in the 
inset. Discrete sources have been excluded from the statistics. 
The solid line represents the best fit of the magnetic field model 
described in the text. The dashed and dot-dashed lines are refer- 
ence models, see text. 

The brightness profile is very sensitive to the value of tj. 
To show this, we fixed the central brightness Iq and we also 
traced the models corresponding to rj = 2/3 (the magnetic field 
is frozen in the thermal gas) and J] - I (magnetic field radi- 
ally decreases according to thermal gas density). The two cases 
are represented in FigfJ] but these two profiles, however, are too 
steep for the data. 

4.2. Numerical analysis of the magnetic field fluctuations 

The analytic model presented in the previous section provides 
a good description of the azimuthally averaged brightness ra- 
dial profile. However, significant deviations of the diffuse emis- 
sion from circular symmetry are observed. These fluctuations 
can be related to local fluctuations of the intracluster magnetic 
field strength and direction around the average. We character- 
ized the power spectrum of these fluctuations following the 3- 
dimensional numerical approaches proposed by Tribble (1991) 
and Murgia et al. (2004). We simulated a multi-scale random 
magnetic field in a cubical box. We assumed that the power spec- 
trum of the magnetic field fluctuations is a power law with index 
n: 

\B k \ 2 cx k-" (7) 

where the wave number k ranges from k m i„ to k max . The sim- 
ulations begin in Fourier space by extracting the amplitude of 
the magnetic field potential vector, A(k), from a Rayleigh dis- 
tribution whose standard deviation varies with the wave number 
accordingto |A^| 2 oc k~"~ 2 . The phase of the potential vector fluc- 
tuations are taken to be completely random. The magnetic field 



is formed in Fourier space via the cross product B(k) = ikxA(k). 
This ensures that the magnetic field is effectively divergence 
free. We then performed a 3D fast Fourier transform (FFT) in- 
version to produce the magnetic field in the real space domain. 
The field is then Gaussian and isotropic, in the sense that there 
is no privileged direction in space for the magnetic field fluc- 
tuations. The power spectrum normalization is set such that the 
average magnetic field strength scales as a function of the ther- 
mal gas density according to Eq.[3] To reduce the computational 
burden, this operation is performed in the real space domain by 
multiplying the magnetic field by Eq.[3] However, to preserve 
the null divergence of the field exactly, this operation should be 
better performed as a convolution in the Fourier space domain 
before the cross product is formed. Moreover, the convolution 
with the radial profile in Eq.[3] alters the input power spectrum 
shape at the edges. Nevertheless, in the specific case of A665, 
the size of the cluster core radius is comparable to the size of 
the computational grid so that the above effects are acceptably 
negligible. 

This 3 -dimensional random magnetic field model is "illu- 
minated" with the relativistic electron distribution in Eq[TJ to 
produce total intensity and polarization synthetic images of the 
radio halo. The magnetic field fluctuates in strength and direc- 
tion from pixel to pixel in the grid. The integration of the total 
intensity straightforward: the contributions of the pixels in the 
computational grid are added together by considering the an- 
gle between the local magnetic field direction and the line-of- 
sight. Indeed, the total intensity at a given direction is propor- 
tional to the integral along the line-of-sight of the magnetic field 
projected on the plane of the sky, B ± . 

The integration of the polarized intensity is not a simply 
scalar sum because of two distinct effects. One is that the mag- 
netic field has random orientations, and the second is that the 
radio wave polarization plane is subject to the Faraday rotation 
as it traverses the magnetized intracluster medium. Therefore, 
the integration of the polarized intensity is performed as a vecto- 
rial sum in which the intrinsic polarization angle of the radiation 
coming in from the pixels located at a depth L is rotated by an 
amount 

A4 1 =Mx|^| 2 , (8) 
where the rotation measure RM is given by 

n e [cm-i]B\\ fyGidl. (9) 



Here By is the magnetic field along the line-of-sight. The combi- 
nation of the aforementioned effects lead to the so-called internal 
depolarization of the radio signal. Finally, to be compared with 
the data, the synthetic images are convolved with a Gaussian 
beam as are the observed ones. The convolution with the beam 
prevents observing the small-scale fluctuations of the radio halo 
and causes a further suppression of the polarized intensity (beam 
depolarization). 

An example of synthetic A665 radio halo is shown in Fig. [3] 
A turbulent magnetic field is simulated in a grid of 1024 3 pix- 
els with a cell size of 1 kpc/pixel. The magnetic field has a 
Kolmogorov power spectrum sloped n = 11/3 and fluctuates 

4 Throughout this paper the power spectra are expressed as vectorial 
forms in £-space. The one-dimensional forms can be obtained by mul- 
tiplying by Ank 2 and 2nk the three and two-dimensional power spectra 
respectively. 
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in the range of spatial scales from A m ,„=4 kpc and A,„ flA =512 
kpc (where A = 2n/k). The central magnetic field strength is 
(Bq) = 1.3 /jG and r\ - 0.47. The simulated magnetic field is pe- 
riodic at the grid boundaries, so the computational grid has been 
replicated to reproduce a field of view of 2048 2 kpc 2 around the 
galaxy cluster center. Left and right panels refer to the total and 
polarized intensity, respectively. Top panels refer to the simula- 
tions at full resolution and show all the fine structure of the ra- 
dio halo. The expected theoretical fractional polarization for this 
power spectrum is about 24%. Middle panels show the images 
convolved at 15" resolution. The radio halo appears smoother 
and the fractional polarization is reduced to about 7%. Finally, 
bottom panels show the synthetic radio images with the same 
noise level as in the observations. In total intensity, only the 
brightest central region of the radio halo are visible, while in 
polarized intensity the halo emission falls below the noise level. 
This example illustrates how radio halos can be effectively po- 
larized, but because of their faintness, detecting this polarized 
signal is a very hard task with the current radio interferometers. 

5. Simulations results 

To constrain the power spectrum of the magnetic field fluctu- 
ations in A665, we compared the synthetic radio halo images 
obtained through the numerical procedure described in Sect. [4] 
and the real radio images and polarization limits presented in 
Sect.[3] Overall, the 3-dimensional magnetic field and relativistic 
electron models depend on the eight parameters listed in Table 
|2] For a power-law magnetic field power spectrum with a radial 
scaling, we have five parameters: (Bo), rj, A mm , A max , and n. As 
pointed out in Murgia et al. (2004), there are a number of degen- 
eracies between these parameters. In particular, the most relevant 
to us is the degeneracies between (Bo) and 77, along with the de- 
generacy between A max and n. Different combinations of these 
parameters lead to similar radio halo total intensity and polar- 
ization distributions. We decided to concentrate our analysis on 
(Bo) (by fixing 77 = 0.47 on the basis of the azimuthally aver- 
aged radial profile presented in Sect. l4.U and on A max (by fixing 
n = 1 1 /3 on the basis of the Kolmogorov theory for a turbu- 
lent medium). Moreover, we fixed A mi „=4kpc, the minimum al- 
lowed by our computational grid. However, we note that a higher 
A„„„ has a negligible impact on the simulation results since, for a 
Kolmogorov spectral index, most of the magnetic field power is 
on larger scales. The random magnetic field is illuminated with a 
population of relativistic electrons with Lorentz factors between 
Jmin - 100 and y max = +00. The electrons spectrum is a power 
law with index 6 = 3, according to the observed radio-halo spec- 
tral index aht - 1 (see Sect.|2]l, and its energy density is in local 
equipartition with the magnetic field energy density. 

Because of the random nature of the intracluster magnetic 
field, the comparison between synthetic and observed images re- 
quires several simulation runs with different seeds for each given 
couple of the fitting parameters (Bo) and A max , this because the 
same set of the magnetic field power spectrum parameters re- 
sults in a different realization of the same synthetic halo. Indeed, 
we chose to simulate 10 equally spaced values of the minimum 
wave number from k m t n = 0.0061 kpc 1 to k m i„ = 0.061 kpc" 1 . 
Correspondingly, the outer scale of the magnetic field fluctua- 
tions ranges from A max - 1024 kpc (the maximum allowed by 
our computational grid) down to A max = 102.4 kpc (i.e. very 
close to the linear resolution of the radio images). For each 
value of Amax, we realized 10 more simulation runs, with dif- 
ferent seeds, in which we determined the value of (Bo) by fit- 
ting the synthetic radio halo brightness profile to the observed 



one. Indeed, we realized a total of 100 simulations by running 
the FARADAY code (Murgia et al. 2004) in the Cybersar-OAC 
computer cluster. 

In Fig.|4]we show the observed radio halo in total emission 
at 25" resolution (panel a) along with the synthetic radio-halo 
images corresponding to three different values A max (panels b, c, 
and d). The synthetic images have the same resolution and noise 
level as the observed one and were obtained with the procedure 
described in Sect. 14.21 First glance, the synthetic images look 
quite similar to the observed one in terms of total extension and 
average brightness. In fact, the top panel of Fig. [4£ shows that 
all three simulations fit the total intensity radial profile remark- 
ably well at 1 .4 GHz of the observed A665 radio-halo image at 
25" resolution. The different lines in the bottom panel of Fig.|4^ 
represent the simulated fractional polarization at different radii 
from the cluster center as expected at 15"resolution before the 
noise is applied. The model fractional polarization is consistent 
with the upper limits derived on the base of the VLA C config- 
uration observation presented in this work. The best-fit central 
magnetic field strength for the three simulations is around (Bo) 
= 1.3 /uG, in good agreement with the value found on the basis of 
the 2-dimensional analytical fit in Sect. 14.11 

Although the examples of the synthetic halos shown in Fig.|4] 
visually appear similar to the observed halo, it is clear that in- 
creases as A, mlv , and the magnetic field power approaches as 
larger scales, the simulated radio halos change shape. In par- 
ticular, when the outer scale of the magnetic field fluctuations 
is close to the observing beam, the halo is smooth and rounded 
(Fig.|4}l). Increasing A max results in a much distorted radio halo 
morphology and in a significant offset of the radio halo peak 
from the cluster centered (particularly evident in Fig.|4j)). The 
simulated clusters are centered on the observed X-ray peak. 

We indeed tried to evaluate quantitatively the value of A max , 
which best reproduce the observed intensity fluctuations of the 
radio halo, as originally proposed time ago by Tribble (1991). To 
do this, we analyzed the residual images obtained by fitting and 
then subtracting, the analytical model in Eq.|6]to both the simu- 
lated and the observed halo images. Since the analytical model 
has a perfect circular symmetry, it is particular well-suited to 
highlighting the halo fluctuations around the average profile. 

In the top panels of Fig|5]we show residuals images corre- 
sponding to the images shown in Fig. |4] with the total intensity 
iso-contours overlapped. The higher A max , the higher the resid- 
ual levels and the larger the fluctuation patches. 

In Fig. [5^ the root-mean-square of the residual images, evalu- 
ated inside the inner ~ 200 kpc from the cluster center, is plotted 
against A max . The solid line and the shaded region represent the 
mean and the standard deviation of the mean derived from the 
simulations, respectively. First of all, the residuals approach the 
noise level as A max approaches the observing beam. This result 
implies that, if A max is smaller than < 100 kpc, the expected rms 
of the halo intensity fluctuations is already below the noise level 
of the observations 07 = 0.06 /xFy/arcsec 2 . As the maximum 
scale of the intracluster magnetic field power spectrum increases, 
the rms level of the radio halo fluctuations increases, reaching a 
value as high as 0.15 /iJy/arcsec 2 for A max = 1024 kpc. 

By comparing the observed fluctuation rms level of 
0.11 /jjy/arcsec 2 with the simulated trend we estimate that the 
outer scale of the magnetic field fluctuations is A max ^ 400 4- 500 
kpc (i.e. k min = 0.014H-0.012kpc~'). In Fig.|5f the offset between 
the X-ray and radio peak is plotted against A max . This indicator 
has a larger dispersion but is much more direct since it does not 
involve a specific handling of the observed and simulated im- 
ages. As observed in Fig. [4] the offset of the radio peak from 
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Parameter 


Value 


Description 




free 


J J It l il\ H_ 1 lltlU > 1 1 L 1 J L; 1 1 1 ul UIL t 1 LIMvl 


n 


0.47 


magnetic field radial decrease, (B(r)) = (B ) (jf^)' 
magnetic field power spectrum index, \B k \ 2 cc kr" 
minimum scale of fluctuation, A mi „ = 2n/k ma t 
maximum scale of fluctuation, A„ ra v = 2n/k mj „ 


n 

^■max 


11/3 
4kpc 

free 


ftnin 


100 


minimum relativistic electron Lorentz factor 


/max 


+oo 


maximum relativistic electron Lorentz factor 


s 


3 


power-law index of the energy spectrum of the relativistic electrons 



the cluster center (here assumed to be coincident with the X-ray 
peak) increases with K max . By comparing the observed offset of 
135 kpc with the simulated trend, it turns out that K max should 
be about 500 4- 600 kpc in agreement with the residual analysis. 

To summarize, the results of the 3-dimensional simulations 
indicate that the radio halo emission in A665 is consistent with 
a central magnetic field strength of about 1.3 fiG. To explain 
the azimuthally averaged radio brightness profile, the magnetic 
field energy density should decrease following the thermal gas 
density, leading to an averaged magnetic field strength over the 
central 1 Mpc 3 of about 0.75 fiG. From the observed brightness 
fluctuations of the radio halo, we infer that the outer scale of the 
magnetic field power spectrum is ~450 kpc, the corresponding 
magnetic field auto-correlation lengtl0 is A# ^100 kpc. 

6. The depolarization of discrete radio sources 

In this section we investigate the effects of the intracluster mag- 
netic field in A665 on the polarization properties of the discrete 
radio sources belonging to the cluster of galaxy itself or situated 
in background to it (see TableO. In particular, it is expected that 
the radio sources whose emission crosses the central region of 
the cluster suffer a higher Faraday rotation, hence a higher beam 
depolarization, with respect to the sources at larger impact pa- 
rameters (see e.g. Murgia et al. 2004). In the left panel of Fig.|6l 
we show a field of view of about 1.5° x 1.5° around A665 taken 
from a 15" resolution image. All the discrete sources brighter 
than 500 juJy/beam (20O-/) are labeled. This severe cut guaran- 
tees a maximum error on the fraction polarization of 0.03. 

For each of these sources we plotted the fractional polariza- 
tion against their distance from the cluster center (right panel 
of Fig|6]l. At large impact parameters, r ± > 2000 kpc, the 
observed radio sources fractional polarization at 1.4 GHz and 
15" resolution oscillates around 0.1. The plot also shows a hint 
that the discrete source seen in projection close to the cluster 
center is more depolarized. 

The continuous line is the expected depolarization trend cal- 
culated on the basis of the best magnetic field model derived 
from the fit of the radio halo image (see Sect|5]l. The dashed lines 
represent the dispersion of the model. The depolarization model 
assumes that all the considered sources are situated at the cluster 
mid-plane and that the intrinsic polarization is fully ordered. As 
intrinsic value we assumed the mean of the fractional polariza- 
tion of the radio galaxies located at more than 2 Mpc from the 
cluster center. 

5 Following EnBlin & Vogt (2003), we define the magnetic field auto- 
correlation length as 

r"|Bii 2 Jfc<a 

As = *-7^ 7^—- (10) 

Jo W k z dk 



Although the relatively large scatter of the data and the sim- 
plicity of the model prevent us from deriving other useful con- 
straints to the magnetic field properties in A665, we note that 
the model prediction broadly agrees with the observed depolar- 
ization trend. Studies at possibly higher resolution and sensitiv- 
ity are required to measure the fractional polarization for many 
more sources so that this method can be applied successfully. 

7. Conclusions 

In this work we presented a new deep VLA observation at 1 .4 
GHz of the cluster of galaxies A665. By combining this obser- 
vation with a previous VLA observation at lower resolution, we 
studied the intracluster magnetic-field power spectrum by ana- 
lyzing the radio halo brightness fluctuations, following both an 
original idea by Tribble (1991) and the numerical approach pro- 
posed by Murgia et al. (2004). 

Our findings are summarized as follows. We simulated 
Gaussian random three-dimensional turbulent magnetic-field 
models in order to reproduce the observed radio halo emission. 
By comparing observed and synthetic radio halo images, we 
constrained the strength and structure of the intracluster mag- 
netic field. We assumed that the magnetic field power spectrum 
is a power law with a Kolmogorov index, and we imposed a local 
equipartition of energy density between relativistic particles and 
the magnetic field. Under these assumptions, we find that the 
radio halo emission in A665 is consistent with a central mag- 
netic field strength of about 1.3 yuG. To explain the azimuthally 
averaged radio brightness profile, the magnetic field energy den- 
sity should decrease following the thermal gas density, leading 
to an averaged magnetic field strength over the central 1 Mpc 3 
of about 0.75 yuG. From the observed brightness fluctuations of 
the radio halo, we inferred that the outer scale of the magnetic 
field power spectrum is ~ 450 kpc, the corresponding magnetic 
field auto-correlation length is 100 kpc. We also find a hint that 
the discrete source seen in projection close to the cluster center 
is more depolarized. The best-fit magnetic model broadly agrees 
with the observed depolarization trend. 

More tightening constraints could be potentially obtained by 
detecting the radio halo polarization fluctuations, not just total 
intensity fluctuations. In fact, the ratio of two former quanti- 
ties, i.e. the fractional polarization, is a very robust indicator 
of the intracluster magnetic -field power spectrum, because it 
only marginally depends on the shape of the energy spectrum of 
the synchrotron electrons and on the equipartition assumption. 
Therefore, it would be very important to improve the sensitivity 
of the future observations in order to detect polarized signal in 
the most radio halo possible. This is a science case in the new 
generation of instruments in radio astronomy. 

Acknowledgements. We thank the referee for very helpful comments that led 
to improvements in this work. This work is part of the "Cybersar" Project, 
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Table 3. A665 radio galaxies used to study signal depolarization through intracluster medium. 



Source 


RA 


DEC 


r± 


FPOL 


S 1.40Hz 




(J2000) 


(J2000) 


(kpc) 


(%) 


(mjy) 


A 


08:33:12 


65:50:41 


2474 


3.0 ±0.1 


222± 7 


B 


08:30:53 


65:36:60 


2429 


2.7 ± 0.2 


122± 4 


C 


08:31:31 


65:52:36 


712 


<0.6 


35± 1 


D 


08:29:51 


65:51:53 


1293 


1.2 ±0.2 


41± 1 


E 


08:30:25 


65:50:42 


639 


< 1.0 


12.8± 0.4 


F 


08:32:02 


65:50:38 


1166 


5.8 ±0.9 


6.4± 0.2 


G 


08:31:40 


65:52:33 


855 


< 3 


3.8± 0.1 


H 


08:28:34 


66:02:17 


3449 


13.8 ± 1.0 


27.2± 0.8 


I 


08:32:17 


66:00:50 


2382 


3.4 ± 0.9 


12.0± 0.4 


J 


08:30:42 


65:58:09 


1456 


< 3 


2.8± 0.09 


K 


08:26:57 


65:42:31 


4745 


6.1 ± 1.8 


19.9+ 0.6 


L 


08:33:34 


65:37:01 


3779 


4.0±2.2 


5.9± 0.2 


M 


08:30:04 


66:13:06 


4259 


3.2 ± 1.9 


5.7± 0.2 


N 


08:34:13 


65:45:25 


3728 


11.3 ± 2.2 


3.7± 0.1 





08:33:11 


65:58:02 


2810 


7.0 ± 1.6 


2.8± 0.09 


P 


08:37:44 


65:13:31 


10203 


5.2±1.5 


1275±2 


Q 


08:34:14 


66:35:52 


8991 


9.4±2.7 


341±2 


R 


08:30:34 


65:51:54 


550 


< 5 


1.05±0.04 


S 


08:31:52 


65:54:53 


1277 


< 6 


0.77±0.03 


T 


08:32:44 


65:50:15 


1946 


< 5 


1.77±0.06 


U 


08:28:51 


65:54:36 


2512 


< 6 


0.81±0.03 


V 


08:29:28 


65:49:08 


1713 


< 5 


1.19±0.04 


w 


08:32:26 


65:34:28 


3313 


< 6 


3.2±0.1 


X 


08:31:08 


65:35:39 


2676 


< 5 


2.31±0.07 


Y 


08:30:09 


65:34:57 


2952 


8.5±2.2 


1.22±0.04 


Z 


08:31:16 


65:47:30 


607 


<7 


0.42±0.02 



Col. 1: Radio galaxy labels; Col. 2 to Col. 3: Source's position; 



Col. 4: Distance from cluster center; Col. 5: Polarization Percentage ; 
Col. 6: Flux density. 
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Fig. 3. Example of simulated radio halo. Left and right panels refer to the total and polarized intensity, respectively. Top panels: 
images at full resolution. Mid panels: images at 15" resolution. Bottom panels: as for mid panels but with noise added. 
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Fig. 4. Observed (a) and simulated (b, c, d) surface brightness images, with FWHM of 25"x 25". Both the observation and 
simulation first contour levels are drawn at 3cr and the rest are spaced by a factor V2. The simulations correspond to three different 
A mflx : 1024, 512, and 100 kpc. In the bottom panels (e), the radial profiles of the radio brightness and the fractional polarization are 
shown. The polarization upper limits are at lcr level. 
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Fig. 5. Observed (a) and simulated (b, c, d) residuals images, with FWHM of 25"x25". Both the observation and simulation first 
contour levels are drawn at 3<x and the rest are spaced by a factor V2. Panel e) and f) show the residuals rms versus A max and the 
X-ray to radio peak offset versus A max , respectively. 
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Fig. 6. Depolarization of the brightest discrete radio sources as a function of the distance from the cluster center (marked by the 
cross). The upper limits are at the 3<r level. The solid and dashed lines respectively represent the mean and the dispersion of the 
depolarization trend expected on the basis of the magnetic field model that best reproduces the radio halo in A665, see text for 
details. 



